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Abstract 

Starting from integrable su{2) (quasi-)spin Richardson-Gaudin XXZ models we derive several prop¬ 
erties of integrable spin models coupled to a bosonic mode. We focus on the Dicke-Jaynes-Cummings- 
Gaudin models and the two-channel (p -|- ip)-wave pairing Hamiltonian. The pseudo-deformation of the 
underlying su{2) algebra is here introduced as a way to obtain these models in the contraction limit of 
different Richardson-Gaudin models. This allows for the construction of the full set of conserved charges, 
the Bethe Ansatz state, and the resulting Richardson-Gaudin equations. For these models an alternative 
and simpler set of quadratic equations can be found in terms of the eigenvalues of the conserved charges. 
Furthermore, the recently proposed eigenvalue-based determinant expressions for the overlaps and form 
factors of local operators are extended to these models, linking the results previously presented for the 
Dicke-Jaynes-Cummings-Gaudin models with the general results for Richardson-Gaudin XXZ models. 


1 Introduction 

Integrable models can be used to describe a wide range of physical phenomena, where the exact solvability 
allows for a numerical treatment beyond the domain of applicability of perturbative and mean-field treat¬ 
ments. One such class of integrable systems is the class of Richardson-Gaudin (RG) systems [1-3]. These 
systems support as many (non-trivial) conserved operators commuting with the Hamiltonian as there are 
degrees of freedom in the system [4] and the eigenstates are given by a Bethe Ansatz wavefunction. A set 
of coupled non-linear equations, the so-called RG or Bethe Ansatz equations, have to be solved in order 
to determine the variables in the wave function. The number of equations scales linearly with the system 
size, in contrast to the exponential scaling of the Hilbert space when diagonalizing the Hamiltonian matrix. 
Another attractive feature of these models is the availability of numerically efficient expressions for overlaps 
and form factors. Traditionally, these expressions can be obtained by Algebraic Bethe Ansatz methods 
and result in Slavnov [5-7] or Borchardt/Izergin determinants [8-11]. An alternative approach has been 
introduced by Faribault et al. [12,13], closely related to the Bethe ansatz/ordinary differential equation 
correspondence [14,15], and later generalized by us towards the full class of RG models [16]. In this ap¬ 
proach an alternative, well-conditioned, set of equations is solved for the eigenvalues of the constants of 
motion characterizing these integrable models. Subsequently, scalar products and certain form factors can 
be written as determinants of matrices whose entries only depend on these eigenvalues [16,17]. Interestingly, 
these expressions can also be seen as partition functions with domain wall boundary conditions [17,18], for 
which similar results were obtained in the context of spin chains [10, 19]. 

Initially, the RG models where formulated in terms of s?.t(2)-algebras describing spin states or fermion 
pairs, but it is also possible to introduce a bosonic degree of freedom by means of a limiting procedure 
[1,20,21]. Two main classes of systems can be obtained in this way. Originally, the class of Dicke-Jaynes- 
Gummings-Gaudin models was obtained from the trigonometric RG model, which includes the Jaynes- 
Cummings [22], the Tavis-Cummings [23] and the inhomogeneous Dicke model [24], all describing the 
interaction between a (set of) two-level system(s) and a single bosonic electromagnetic mode. Another 


1 


class of systems considers the coupling of an integrable two-channel {p + ip)-w&\e superfluid to a bosonic 
mode, which can be obtained as the limiting case of a hyperbolic RG model [21]. This model was initially 
introduced by Dunning et al. [25] and can be shown to be equivalent to a model which couples Cooper pairs 
to condensed molecular bosons [26] . These two classes were later obtained as two distinct cases of integrable 
Hamiltonians containing a bosonic degree of freedom, starting from a variational approach [27]. 

The derivation of form factors and overlaps in the eigenvalue-based formalism [16-18] depends heavily on 
the existence of a dual state and therefore a (dual) highest weight state in conjunction with a lowest weight. 
This makes the generalization of results for the su(2)-models towards models containing a bosonic degree of 
freedom far from straightforward, because the hw{l) algebra of a bosonic mode is non compact and therefore 
lacks a highest weight. Tschirhart and Faribault recently showed how determinant expressions could be 
found for the form factors of the Dicke-Jaynes-Cummings-Gaudin models by means of the introduction of 
an intricate alternative Algebraic Bethe Ansatz [18]. In the present paper, we obtain these expressions as 
a limiting case of a renormalized pseudo-deformed spin model. The pseudo-deformation scheme [28] was 
originally proposed as a way to shed light on the singularities arising in the RG equations by connecting the 
spin models to purely bosonic models. Here this method is used as a way to obtain a bosonic algebra as the 
contraction limit of a su{2) quasispin algebra. The connection can be made adiabatically in a controlled 
fashion [29] and allows for an extension of the su(2)-based integrable systems towards those containing a 
bosonic degree of freedom. The purpose of the present paper is twofold. First, we show how the eigenvalue- 
based formalism for the Dicke model, introduced by Tshirhart and Faribault, can be generalized to the 
extendend (p + ip) (and related) models in the same way we recently generalized the su{2) XXX results 
to XXZ models. Second, we present the pseudo-deformation scheme as a unifying framework in which 
all properties (eigenstates, eigenvalue-based variables, form factors, etc.) of both the Dicke and extended 
(p + ip) model can be derived in a simple and straightforward way. 

In section 2, the necessary preliminaries for Richardson-Gaudin models will be reviewed, after which 
the pseudo-deformation scheme is discussed (section 3). The connection with the bosonic models is then 
made explicit in sections 4 and 5 and determinant expressions are presented for overlaps and normalizations 
in section 6. It is shown how the results for XXZ RG models can be generalized to these models, further 
extending the description of Richardson-Gaudin integrable models in terms of ’eigenvalue-based’ variables 
and unifying some previously presented results [16, 18]. 


2 Richardson-Gaudin models 


2.1 Definitions 


Richardson-Gaudin models are defined by a set of n mutually commuting conserved charges [30,31] parametrized 
as 


R, = S? 


k^i 


t* 

2 


(Sis. 


+ SkSj) + Z,kS°S° 


( 1 ) 


with the set of operators {Sj, Si, S^} {i = l...n) spanning a set of n independent su{2)i (quasi-)spin 
algebras 

[5°,5^] [Sf,Sk]=-S,kSk, [4,5fe]=2d,fc,5°, (2) 

with irreps \si,p,i) associated with each separate algebra 5^(2)^. These algebras can represent genuine 
spins, general {2si + l)-levels by means of a Schwinger representation, or fermion quasispin pairs in a 
pairing model [32]. The introduction of these operators allows for an algebraic formulation of these models, 
independent of the underlying physical interpretation. 

The constraints on the X- and Z-coefhcients for which these operators commute mutually were obtained 
by Gaudin [1] and result in a set of equations defining a Gaudin algebra as 


Xij — Xji, 


Zij — Zj i , 


XijXj]^ Xii^(^Zij -j- Zj]^) — 0 , 


(3) 


which have to hold Vi j ^ fc = 1,..., n. Multiple classes of solutions for these equations have been found 
[1,30], where each class considers and as odd functions of a set of parameters {ci} = {ci, £ 2 , • ■ •, Cn}? 
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such that the Gaudin algebra is defined as = X{ei,€j) and Z^j = Z{ei,€j). A model is said to be RG 
integrable if it has as many conserved charges commuting with the Hamiltonian as degrees of freedom [30,33] , 
so an integrable Hamiltonian can be constructed as a linear combination of the conserved charges (1), 
commuting with all these operators by construction. 


2.2 Diagonalizing integrable Hamiltonians and eigenvalue-based variables 

Since the conserved charges (1) commute mutually, they share a common set of eigenstates. However, since 
any arbitrary product state can be created by either acting with creation operators on the lowest-weight 
state or by acting with annihilation operators on the highest-weight state, two separate representations 
are possible for each eigenstate. These will be discussed simultaneously here. After extending the Gaudin 
algebra (3) by adding a set of variables (also called rapidities) Xa such that Xi^ = X{ei,Xa), the following 
creation/annihilation operators can be defined as 

n n 

^ X,^sl S^ = Y, X,^S,, (4) 

i—1 i—1 


fully determined by a single (possibly complex) variable Xa, where Xia = X{ei,Xa) extends the Gaudin 
algebra. Eigenstates of the operators Ri can then be constructed by the repeated action of generalized 
creation/annihilation operators on an empty/fully-filled vacuum state 

C N \ /M-N \ 

iV'')=(n1^')- (5) 

where we have defined the particle-vacuum state \9) = jsi,—Si) and the hole-vacuum state \9') = 
G?=i \si,Si). The number of excitations N is restricted to N < M = Both states \tf}) and ['!/'') can 

represent the same eigenstate of Ri, in which case they are referred to as dual states. These states can be 
shown to be eigenstates if the RG equations 

n N 

1 + g'^SiZ.a - g'^ Zpa = 0 (particles), 

2—1 

n M-N 

-l- g'^SiZ^a' - g ^ Zp'a'=Q (holes), (6) 

2—1 

are satisfied. The resulting eigenvalues of the conserved charges Ri are then given by 


N 


— di I 1 g ^ ^ ^ict T g ^ ^ ^ikdk 

k^i 


Ct — l 

M-N 


r' = di I 1 - g y] Zia' + 5 T Zikdk 

k^i 


(particles). 


(holes). 


(7) 


Ot' — l 


Following on recent work on the XXX RG model [12,13], we showed [16] how an alternative description 
of the Bethe Ansatz states could be obtained starting from a new set of variables 


N M-N 

K, = ^Z,a, A'= y] Zia', Vz=l...n. (8) 

ct—1 ol' — 1 

It can be seen that each eigenvalue (7) is only explicitly dependent on the parameter A^ (particle) or A' 
(hole). This has led to the denomination ’eigenvalue-based variables’. Since the eigenvalue of an eigenstate is 
independent of its particular representation, a correspondence between A^ and A' can be found by equating 
the eigenvalues = r(, leading to 

gK = dAj + 2, Vz = 1... n. (9) 


3 


So, if the eigenvalue-based variables are known in one representation, those for the dual representation 
immediately follow. Remarkably, the overlaps with non-interacting states (such as single Slater determinants 
for the quasispin picture), as well as the normalization and several form factors can also be determined solely 
from these variables (see the next subsection). 


2.3 Determinant expressions 

When performing calculations with Bethe Ansatz states, it is often convenient to expand them in a basis 
set. The wavefunction can be expanded in the complete set of basis states 


N' 

iw})=n( 5 !) ‘| 0 ), 


( 10 ) 


with the integers [Ni,... ,Nn] = [{IVi}] a partitioning of the number of excitations N over the number of 
levels n. The expansion is given by 


N 


n 


. . Ni 

E n (-^0 1 ^)’ 

nm] ■ ^ 


OL — l \ 2=1 


2 = 1 


with expansion coefficients given by the permanent of a matrix [34,35] 


and 


r^N 


N, 


Nr. 


( 11 ) 


( 12 ) 



• ^iiOLx 


Xi^ai ^ 



. ^iiOL2 

Xijyory 

Xijyay 

(13) 

^ixOiN 

^ixOLN 

XirydM 

Xi^ay, / 



For donbly-degenerate models (si = 1/2, Vi), the occupation W of each level is either 0 or 1 and this 
permanent can be rewritten as a determinant [16], which can efficiently be evaluated numerically. Labelling 
the set of occupied levels (W = 1) as {ia} = {h, ■ ■ ■, *Ar}, the overlap of a Bethe Ansatz state with 


|{*a}) = Aat}) 


N 




V,a=l 


(14) 


is given by 


with 


= ("(*1 ■ ■ • w }|{ a : a }) 



det J 


Jab — 


- E 


N 

c^a 




if a = 6 
if a b 


(15) 


(16) 


Here the Gaudin algebra (3) was extended by means of an arbitrary gauge variable as introduced in [16]. 
It is worth noting that the overlaps ({H, • • ■ )*Af}|{a^a}) (Eq. (15)) are implicitly independent of the gauge 
variable e^, so the freedom left in the choice of this variable can be exploited to obtain the most simple 
expression possible for each realization of the Gaudin algebra, usually by setting = 0 or oo. This 
expression for the overlap, combined with the existence of the dual representation, are the key ingredients 
for the following results. The overlap between a state and its dual state is given by 




{e'\ 




(17) 
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which can also be seen as the overlap of a Bethe Ansatz state with M excitations defined by the set of 
variables {Xa} U {xa'} with the hole-vacuum 


M 

(K}iw)=(0'i n 


(18) 


For this state, the eigenvalue-based variables appearing in the diagonal elements of the overlap matrix are 
given by 

N M-N 2 

^ia + Zia' = Ai -|- A( = 2Ai H—, (19) 

a—1 a' — l ^ 


where we have used the relation (9) between these variables for the dual representation. The overlap (18) 
can then be written as 


({^'a}|{^4) 


rig ^ra 


det J, 


( 20 ) 


with 


2Ai + I ~ Zik + if j — j 

Xij if * j 


( 21 ) 


Note that this does not depend explicitly on the dual state, but only on the original state through the terms 
Ai in the diagonal elements. Once this overlap is known, together with the overlap of both representations 
with an arbitrary basis state, all information is present to calculate the normalization of both states, as 
originally shown for the RG XXX model [17]. 


3 Pseudo-deformation of the quasispin 

The results reviewed in the previous section can be generalized to integrable systems containing a bosonic 
degree of freedom by replacing one of the su(2) spin algebras with a bosonic hw{l) algebra. There are 
multiple approaches possible for this process of bosonization, of which the Holstein-Primakoff transformation 
[36] is arguably the most well-spread. Here we have opted to use the recently proposed pseudo-deformation 
scheme [28] as a way to obtain bosonic commutation relations because this method provides an adiabatic, 
and therefore controlled, mapping of the hard-core bosonic su{2) algebra on a genuinely bosonic algebra. 

A pseudo-deformed su{2)^ algebra can be defined as 

[5°(C),5f(0]=^^(0, [S°{0,S{0]=-S{0. [SH0,S{0] = 2{iS°{0 + i^-l)s) (22) 

with ^ S [0,1] the pseudo-deformation parameter and s the original (^ = 1) su{2) irrep label. This definition 
can be interpreted as providing a linear interpolation between two known limits: ^ = 1 gives rise to the 
original su{2) quasispin algebra, while ^ = 0 results in a (unnormalized) bosonic hw{l) algebra. This latter 
limit was also termed the contraction limit of the algebra [37]. The nomenclature pseitdo-deformation was 
originally proposed because this algebra can be reduced to a canonical su{2)^ algebra 

[^°(0,^^(C)] = ^^(0, [^°(0,^(0] = -^(0, [^ne),^(C)] = 2A°(0, (23) 

by defining 

= A{0 = ^Si0, A°(e)=^°(O+(l-0A (24) 

except for the contraction limit ^ = 0. In this limit, the following operators 

= b^b = S°{0) + s, (25) 

close the bosonic hw{l) algebra 

[b'^b,b^] = b\ [b'^b,b] = —b, [6,5^] = !. (26) 
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The irreducible representations of the A(^), algebra are labeled by s(^) = s/^. The inter¬ 

pretation behind this is a gradual increase of the effective multiplicity (2s(^) -I- 1) of the su{2)^ irrep with 
decreasing It should be noted that only discrete values of ^ (with n = 2s, 2s -|- 1,...) give rise 
to unitary irreps. Nevertheless, this is not problematic because the theory of RG integrability is not based 
on matrix representations (with integer dimensions), so the parameter ^ can be regarded as a continuous 
variable. 

This construction has led to a numerical solution method for the RG equations ( 6 ) by solving the 
equations adiabatically from the contraction limit (^ = 0) to the ^ = 1 case [28, 38]. Because the 
{A°(^), Al(^), y4(^)} operators span a canonical su{2)^ algebra, the set 




k^i 


+ Ak{i)A^ (^)) -I- ZikA^ (^)A[ 


.(0 


Vz=l,...,n, (27) 


remains in involution for every value of ^ 7 ^ 0 provided the Gaudin equations (3) are satisfied. The Bethe 
Ansatz state in the particle representation for these operators is then given by 


N / n 

i^)= n 

if the pseudo-deformed RG equations 

n N 

l + = 0 , Va = l...Af, (29) 

2 — 1 



are satisfied. The key observation here is that the set of coupled non-linear RG equations ( 6 ) reduce to 
a single uncoupled equation in the contraction limit ^ = 0 , formally equivalent to the secular equation of 
the pp-TDA [39]. This equation can be straightforwardly solved numerically, after which these solutions 
can be adiabatically connected to the RG equations of the original problem ( 6 ) by slowly tuning ^ 1, as 

proposed and illustrated in [28,38]. 


4 Integrable models containing a bosonic degree of freedom 

In the previous section a fully-bosonic integrable model was obtained by deforming all su(2) algebras 
simultaneously into hw(l). A natural digression would be to consider a situation where only one of the spin 
algebras is pseudo-deformed, leading to an interacting boson in the contraction limit. Since the bosonic 
limit is a singular limit of the algebra, care has to be taken to obtain finite results in the contraction 
limit. However, if a model is found for which no singularities arise in this limit, an RG integrable model 
is obtained. Several integrable models are known which also contain an interacting boson, of which the 
Dicke-Jaynes-Gummings-Gaudin model (henceforth referred to as the Dicke model), is the best-known. We 
will detail the derivation for the Dicke model, and summarize the key results for the extended (p + ip) model, 
due to the large similarity with the Dicke model. 

4.1 The Dicke model 

The Dicke Hamiltonian [24] is given by 


m m 

H = eob^b + ^ e,S° + 5^ (^1^ + *^*^0 
2—1 2—1 

and describes a set of m two-level systems (si = 1 / 2 ) interacting with a single mode of the bosonic field, 
represented by a photon with energy eg. The connection between this model and the XXZ RG systems was 
already made by Gaudin [1] and later extended by Dukelsky et al. [20]. The derivation presented here is 
similar to the one by Dukelsky et al., but differs in our choice of bosonization scheme [30]. Starting from 
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the constants of motion for a set of m + 1 spin systems, we label them i = 0,1 ,... ,to and exchange the 
su{2) algebra labeled i = 0 by a pseudo-deformed algebra su{2)^. Starting from the trigonometric Gaudin 
algebra [31] 


Xij — 


1 + »??)(1 + ??|) 

Vi - Vj 


1 + r]^r]j 

Vi - Vj 


(31) 


and taking the limit t oo, a Gaudin algebra is obtained determined by 


Xok = a/1 


vl 


Zok = Vk, 


Xik — 


^ (1 + Vi)i^ + v]) 


Zih — 


1 + ViVj 

V^ - Vj Vi - Vj ' 

The conserved charge associated with the deformed algebra is given by 


(32) 

(33) 


Rq{0 — 

k^O 

m 

= A°(c) + ff^ 

A: 7^0 


^Xok {A^OSk + AiOSl) + ZokA°iOS? 
\\S^k {A\OSk + ^(04) VkA\£,)S^ 


(34) 


It is now possible to define a ^-dependent coupling constant as ^ and rescale the variables 

T]k = ^63-^ 111® contraction limit (^ « 0) the Gaudin algebra now reduces to 


which is related to the parametrization proposed by Dukelsky et al. [20] . By making use of this parametriza- 
tion, the Dicke Hamiltonian can be obtained (up to a diverging constant) as 

m m 

< = eoRoi^ ^ 0) = eob^b + ^ +gJ2{sjb + S,b^) (36) 

i=l i=l 


in the contraction limit. A similar procedure leads to 

1. the other constants of motion Rf [30] 


R 


D 


= eoR^{^ ^ 0) = (eo - e,)5° - 0(^6 + S,b^) -2G^Y.— 


k^i 




^{SjSk- 


• s.sl) + s°s° 


(37) 


2. the Bethe Ansatz for the Dicke model [40] 

where the prefactor can be absorbed in the normalization, and 

3. the RG equations [40] 


m N ^ 

(eo-cca)-2G^y +2G^y-= 0, Va = l...A^. 


(38) 


(39) 
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Remarkably, it is also possible to transform the equations for the rapidities to an equivalent set of equations 
for the set of eigenvalue-based variables provided Si = 1/2,Vi = [12,41]. For the Dicke model, 

these variables are given by 


N 


A, 


= E 

q ;=1 


= 1 , 


. m. 


(40) 


These determine the eigenvalues of the constants of motion (37) as 


Ri 


D I 


N 


e,-eo) + 2G^Y. -E- 




k+i 


£i ~ Co) + 20"^Ki — ^ ] - 




\i’) 


and satisfy the coupled quadratic equations 


rri . . 

= Ar_A,(e,-eo) + G2y Vz = l...n. (41) 

These equations can either be determined starting from the RG equations for the Dicke model [12,18,41], 
or by taking the contraction limit of the eigenvalue-based equations for the trigonometric RG model (see 
Appendix A). The quadratic equations can easily be solved numerically, whereas the RG equations become 
singular near the so-called singular points [42,43]. Similar to the SM(2)-based RG models, it is possible to 
express the overlap and form factors of the Dicke model in the eigenvalue-based variables only, circumventing 
the need to calculate the singularity-prone rapidities [18] (see also section 6). 


4.2 The {p + ip)-'wacve pairing model 

A similar procedure can be used to prove the RG integrability of the {p + ip)-wave pairing model coupled 
to a bosonic degree of freedom. This model was introduced by Dunning et al. [25] as an extension of the 
integrable fermionic p^ + jp^-pairing model [44-46] . Lerma et al. consequently showed how this model is 
given by the limit of a hyperbolic RG model [21], which we will reformulate by making use of the pseudo¬ 
deformation. In addition, we will show how the model fits within the eigenvalue-based language. The 
Hamiltonian is given by 

H = 6b'^b + ^^clck-^{k^-iky){k'^+iky) y cJ^cl^Ck'C.k' 

k k^±k' 

“ y E )clcclk^ + ^-c-) > (42) 

k 

and was shown to be integrable if 

5 = -F^G, K = FG. (43) 

By making use of the quasispin formalism [32] and absorbing a phase in the quasispin operators, the 
Hamiltonian can be rewritten as 

m m _ m 

H = Sb^b + Y,^kS'i-G Y. ^fke^^SlSk'-KYV^{slb + Skb^y (44) 

k—1 k,k' — l k—1 


which can again be related to a SM(2)-based RG model. Starting from a hyperbolic Gaudin algebra [31] 


^ _ o ^ _ ei + ek 

uXij — Z , Z/ij — . 




^k 


( 45 ) 














for (m + 1) levels and exchanging a single su{2) algebra with a pseudo-deformed A(^)-algebra, the Hamil¬ 
tonian can be obtained as a linear combination of the conserved charges in the contraction limit. Labeling 
this single pseudo-deformed algebra as ’0’ and renormalizing the coupling constant and bosonic energy level 
eo as 


So 


C 2 

eo = 

2 so 


(46) 


the constants of motion become 


m 



= lim ^ 


^- I 

lim ‘fRM) =Y. ^ -5*5^ + S.Sl) + (l±(isfS2 

C-*-o ? “ — Efc \ / e* — fife 

k^i 


?7o_ 

r^i 




ei 


(47) 


(48) 


with N the number operator counting the number of excitations. Note that i?i(^ —)■ 0) = 0, but 
remains finite for the whole range of ^ and results in a non-zero conserved operator in the contraction limit. 
These operators are the building blocks for the Hamiltonian (44), similar to the results presented for the 
fermionic Px + ipy pairing model [45] . We can take the linear combination 


m m 


I K + b^b + y/^k^k' 

k'^k / k^lk'^k 






^0 


m m 


5! 


k 5 


/c=l 




(49) 


then add the Casimir operators for each algebra su(2)k times Ck, and finally introduce the number operator 
(47). The resulting Hamiltonian now becomes exactly Eq. (44) 

m m 

H = r^lb^b + Y,{^ + N)ekS°+ ^ 

k—1 k,k' — l 

m 

+ ^0 ^ + Skb^^ - rj^N. (50) 


The integrability condition (43) arises naturally from the parametrization of the Gaudin algebra. The 
Hamiltonian studied in [26] can be obtained by taking the linear combination 


mm m 2 ^ 

^ = E < = E ^ ■ E + E +^)- 


m 9 

— Vo 

i=l i=l V i-i 

The Bethe Ansatz states for these models can also be found from the contraction limit as 

N / m ,— ' 

^k ^OL ^0 




Q;=l 


k^\ 


with resulting RG equations 


N 


-E 


1/3 -I- Xo 


2 

Vo 1 + Xo 

K -^ - 

r,, ^/3 

k—1 /3#a 


= 0, Va = 1... TV. 


(51) 


(52) 


(53) 
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Using the techniques from [16] (see Appendix A), these can again be shown to be equivalent to a set of 
quadratic equations in the variables 


N 


Ai = y -, Vz = 1... m, 

a—l 
N 9 

A. = E 


(54) 

(55) 


which have to satisfy 

= —N(m — N) — 2KKi + 2ril (aq + ^ ^3 — Aj), Vz = 1... m, (56) 

m 

2Ao = ^ A, + 2kN. (57) 

i=l 

The name ’eigenvalue-based variables’ is apt since a single variable A^ fully determines the eigenvalue of a 
single constant of motion R^. If the rapidities satisfy the Richardson-Gaudin equations, we obtain 


<IV’) 


1 

2 



Q!=l 


“h ^oi 


Ci XcK 


2 

Vo 


ei 


m 


+E 


Cj + Efc 
^2 


IV’) 


1 

2 



m 


+E 


Ci + efc 
Ci ~ Cfc 


IV-)- 


(58) 

(59) 


5 Investigation of the adiabatic connection 


Since the pseudo-deformation provides an adiabatic connection between the bosonic algebra and the qua¬ 
sispin algebra, there is a similar connection between the Dicke model and an XXZ RG model. The adiabatic 
connection can be investigated in order to shed some light on the connection between the two models. When 
deforming a single level, the RG equations for arbitrary ^ are given by 

m N 

l + gso(C^oa + |E^“-5E^/3a=0, Vo = 1 ... TV. (60) 

2=1 0i^O' 


Inserting the Gaudin algebra and the parametrization for the Dicke model in these equations, the RG 
equations at arbitrary values of ^ are given by 


{eo- 


2G^ + ^CiEc/so 
ii- Ea 


N 

E 

/3#« 


2G2 + ^EpE^/so 


Efi — Er, 


= 0 , 


Va = I... V. 


(61) 


For ^ = 0, these equations reduce to the equations for the Dicke model, while for ^ = 1 we obtain RG 
equations for a trigonometric spin model 


{eo- Ea)-J2- 


2G^ -f CiEct/ sq 
Ci - E„ 


N 


2 G^ -|- EpEa/so _ Q 


/3#a 


Efi — E„ 


Va = 1... V. 


(62) 


Interestingly these equations can also be seen to reduce to the regular Dicke model equations for an infinite 
Sq. This is easily understood as it corresponds to an infinite degeneracy and thus a bosonic mode, so the 
pseudo-deformation scheme would not change the physics of the problem. When solving equations (61) for 
arbitrary it can be seen that singular points occur, similar to those occurring in the RG equations for 
^ = 1. In these singular points, multiple rapidities Xa coalesce with a single-particle level e^, leading to 
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diverging contributions to the RG equations. These divergencies cancel mutually exact, but have prevented 
straightforward numerical solutions of the RG equations for a long time [43,47]. The occurrence of singular 
points can be linked to the Pauli exclusion principle [28], and no singular points occur in a fully bosonic 
limit. However, since we are only dealing with a single bosonic mode coupled to spin modes, singular points 
will remain in these equations at every value of As an illustration, several solutions to the equations 
have been given in Figure 1 for all values of the deformation parameter ^ = 0... 1, showing the qualitative 
behaviour of the solutions between both limits. As can be inferred from Figure 1, the qualitative behaviour 
of the connected eigenstates in the Dicke model and XXZ RG model can be quite different. For instance, 
the parameters used in the model are such that the Dicke model is in the weak-coupling limit (small G). 
Nevertheless, the connected XXZ RG state in the second column is clearly a collective state in the Gooper 
pairing regime [4], as indicated by the complex-conjugate rapidities. The qualitative difference is accentuated 
by the presence of multiple singular points along the adiabatic path of the pseudo-deformation. 











0.0 0.2 0.4 0.6 0.8 1.0 


0.04 - 
0.02 - 

0.00 - 

- 0.02 - 
- 0.04 - 

- 0.06 - 1 - 1 - 1 - 1 - 

0.0 0.2 0.4 0.6 0.8 1.0 


Figure 1: Several examples of the connection between the Dicke model and the XXZ model for N = 6 
excitations. Each column depicts the connection between different eigenstate of the models. The parameters 
are chosen such that the Dicke model (left side of each column) is in the weak-interaction regime (G = —0.1, 
eo = 1, {ci} = {2, 3,4,..., 12} and sq = 1/2). Although the rapidities for the Dicke model approximately 
coincide with the single-particle levels, they do not end up in the weak-coupling limit for ^ = 1 in the XXZ 
RG model. 


6 Deriving eigenvalue-based determinant expressions for form fac¬ 
tors and overlaps 

In the previous section, it was shown how many results for integrable systems containing a bosonic degree 
of freedom can be obtained in the contraction limit. However, if we wish to obtain expressions for normal¬ 
izations and form factors starting from the RG models, we face the problem that the bosonic level has no 
state of maximum occupation number. A hole-vacuum now cannot be defined since this would contain an 


11 

























infinite number of bosonic excitations. This problem was first envisioned for the Dicke model by Tschirhart 
and Faribault [18], who devised an alternative formulation of the Algebraic Bethe Ansatz, introducing a 
pseudovacuum, which allowed for a description in terms of the eigenvalue-based variables for the Dicke 
model. As we will show, the pseudo-deformation scheme allows for a simpler derivation of these results by 
means of a renormalization procedure. Starting from the eigenvalue-based determinant expressions for the 
XXZ RG models [16], the dual state can be defined for any Bethe Ansatz state. By dividing the overlap 
between the Bethe Ansatz state and its dual state by the overlap between this dual state and a reference 
state (see section 6 . 2 ), a finite expression is obtained in the contraction limit, corresponding to the results 
previously presented [18]. This approach can then be extended to the (p -|- *p)-wave pairing model. 


6.1 Overlap with non-interacting basis states 

A determinant expression for the overlap of an arbitrary Bethe Ansatz state in the Dicke model with a 
basis state (10) will first be derived. The expressions for these overlaps do not depend on the existence 
of a dual state, nor a hole-vacuum state, so the inclusion of a bosonic level into the existing fermionic 
su{2) models in the pseudo-deformation is well defined. The main difference with the previously-considered 
spin-1/2 models is the occurrence of the bosonic level with arbitrary occupation number. The expansion 
( 11 ) with permanents as expansion coefficients ( 12 ) still holds, but since multiple columns of the permanent 
can be equal, these can not immediately be rewritten as determinants. This problem can be avoided by 
introducing a limiting procedure. The permanent of a coefficient matrix with distinct columns can be 
rewritten as a determinant (15), so we will introduce a matrix with different columns and consider the 
limit where multiple columns become equal. The overlaps can then be found from a two-step limiting 
procedure: first the permanent with multiple equal columns can be reduced to a determinant, after which 
the contraction limit can be taken. 

Assume a model where all spins Sj = 1/2, except for one level (again labeled 0) with an arbitrary large 
degeneracy. The Bethe Ansatz can then be expanded in a set of basis states as 


N / 771 \ ^ /N-No \ 

n xoT^si+j2x,^si |0) = E ( 5 ^) ” n < 1^) > 


[{V,}] 


(63) 


with the vacuum state now containing the lowest-weight state [sq,—sq) for the degenerate level. The 
expansion coefficient is given by 





-^OCKl 


XiM-No<^l \ 

1 

^Oa.2 


^Oa.2 

Y 

Y 

-^iN-NQOL2 







\ -^OCKiV 


-^Oaiv 

^ilOLN 

^IN-NqOLn ) 


No 


(64) 


In the following, it will be shown how this permanent can be written as the determinant of an (N — Nq) x 
{N — A/))-matrix for the trigonometric realization of the Gaudin algebra leading to the Dicke model. For 
this realization, the Gaudin algebra elements (35) associated with the bosonic level are given by 


Xoc. = V^ + ^a 


lim 

r]o^oo 


+ Vo)i^ + vl) 

Vo - Va 


(65) 


Instead of immediately taking the limit ryo —> oo, it is possible to first evaluate the permanent for arbitrary 
ryo and later take this limit. We now have a parametrization with a free parameter tjq and wish to evaluate 



/ ^Oai 

^Oqi 


Y 

^iN-NQOil 

r/o—>-oo iVo! 

^0a2 

-^Oaiv 

■^iia2 

^iN-NQOi2 



■ -^OCKiV 


Y 

^iN-N^aN 


No 


( 66 ) 
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Instead of using the same parameter tjq for each column and taking the limit for each column simultaneously, 
it is possible to introduce a different variable Q (replacing 770 ) for each column i, and taking consecutive 
limits to infinity of these parameters. This reduces the problem to the evaluation of 


= TVT lim 1 ™ per 

iVo! Cl^OO CiVo^oo^ 


( -^Oiai 



s 

0 

1 

^OlOn.2 

^0NQa2 

^ilO'2 


\ ^0l QiV 

^OnqCun 

^iiaN 

Y 

^IN-NqOCN 




(67) 


No 


with 


^OiCi — 


Vii+cm+vi) 


( 68 ) 


Ci-Va 

This is the permanent of a matrix where each matrix element satisfies the Gaudin algebra (3), for which a 
determinant representation exists 


1 




= lim lim - , - 

iVo! Ci^oo Uo -00 ^N-No + 


det J, 


with J defined as 


Jii — 


1 i+mric i+ViVk 

U ^ r}i-n„ n-rii. ‘ — J 


Vi-Vj 


iii^j 


(69) 


(70) 


This is simply the application of Eq. (15) with the auxiliary level tending to infinity. For the problem at 
hand, we can identify four different sectors in the matrix. If we associate the indices i,j with levels tending 
to infinity and the indices a, b with the finite levels, J can be expressed as 


Jij 

J^a 

Jai 

Jab ^ 


J = 


where the matrix elements of these sectors are given by 


{ j^. I l+CiTjc. _ sr^No 1+CiCj _ sr^N-No 1+CiVo. ,■ _ „■ 

Vi -r Z^a=l Ci-»)Q 0-Cj l^a=l Ci-r]a ‘ “ J 

c.-o 


iii^j 


Jab — 


l + VaVa Y^N-Nq l + rja-qt \^N l + 7]aCi 'f _L 

\ la \ 2-^a—l Ci—Va 2-^b^a rja — rji, 2-^i—l rja—Ci 

\ Vn—Vh 


\i a 


Jin. — Jni — 


\/i + CVi + ??g 


(71) 

(72) 

(73) 

(74) 


G - Va 

Absorbing the factors (1 + ~ \fCi from the prefactor into the first A"o columns i and the first Nq rows 

i and taking the subsequent limits to infinity as 


the overlap can be written as 


Cl > C2 > • • • > Cwo: 


= 1 Jii/^idet J. 


(75) 

(76) 


Here all matrix elements were redefined as (with lim denoting the subsequent limits to infinity) 


Jij — 


lim 


lim 




I i+Ciaa _ y^-^o i+CiCfc y^.iv- 

Ci — 7ja 2-^k^i Ci~Cfc 2-^a— 


N—Nq 

1 Ci-Va 


1 (1 + ^ 2 ) 174(1 + . 2 ) 1/4 


C.-C: 


if i = j 

Hi 


(77) 
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^ab — 


^ l + ’?a’?Q ^^N—Nq l-\-r}arjb l + r)aCi 

— rih 2^i—l 77a—Ci 


n -L i-+r]aVc. Y^iV-- 

' ^-^cx—1 rja — Tjoc 

yi + ^aVi^ 


Va-Vb 


Jin — Jai — lim 




\/i + CVi + ^g 

Ci - Va 


if a = 6 

if a 7 ^ 6 

(78) 


(79) 


All limits are straightforward except the limit of the diagonal elements in the first sector, which needs 
a few algebraic manipulations 


Jii = lim 


= lim 


= lim 


= lim 


= lim 


1 


1 

1 

1 

1 


* 1 + CiVa _ f + CiVa _ 1 + QCj 

G - Va Ci - ria Ci - Cj 

N N—Nq Nq 


c.+E-)-- E ’<.-EV^ 


ot—1 a—1 


^ ^ 0 _ 1 -I- _ 1 -I- d d 

c, + E>7.- E i-E^-E^ 


q;=1 a—1 


j<i 


N N-Nq 

q;=1 a—1 j<ii j^i 


Ki + X! ^3 


j>i 


Now all matrix elements can be determined as 


^ Va—Vb ' ’ 


Jij — 


i ii i = j 

0 ifi^j, 

^ia — Jai — 0 - 


(80) 

(81) 

(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 


In summary, the first Ag diagonal elements become 1,2,..., Ag, while the off-diagonal elements in the first 
Nq rows reduce to zero. The expansion coefficient can then be rewritten as 



■ 1 .. 

. 0 

0 

, ^ UaV^ + ^l J . 

Hn.}] ^ bet 

riz^o V1 + Ti 

0 .. 

. A^O 

0 


0 ... 

Jab 


rig \/l +^a 

n*^o 


det[J], 


( 88 ) 


with the {N — Nq) x {N — Nq) matrix Jab defined as in Eq. (85). This holds for arbitrary parameters 
{r]a, a = 1.. .N — A^o} and {? 7 q, a = 1... A^}. In order to obtain expressions for the Dicke model, it is 
possible to again introduce the parametrization (35) and take the contraction limit. This immediately 


results in 





(Nq-, {ia}\ilj) = 

= V^(-G)'^-'^“det[J] 

(89) 

with 





(spN 1 

J j Z-^a — 1 —X 

^ab — A X 

1 

V—vAT—A^o 1 ’f t 

if a 7 ^ 6 

(90) 
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A similar derivation for the (p + ip)-wave pairing Hamiltonian (44) is given in Appendix B, where it is shown 
that the overlap of a Bethe Ansatz state 


N 






Q!=l 




|0>, 


with a non-interacting state (10) is given by 


{No-, {lam = 


v'fVoI det J 


with J defined as 


Jah — 


■ 1 Ea+Xa V^Af-Wo 1 Ei^+Ei^ A^O + l U _ J, 

2 Z^a=l ei„Z^c/a 2 e ,- e , 2 Hfl — 0 




if a ^ & 


(91) 


(92) 


(93) 


6.2 Normalization 

The procedure in the previous section can also be used to obtain the normalization of the Bethe Ansatz 
states for these models. Assuming a dual state |'0O exists for each eigenstate |^), which is identical but has 
a different normalization, we have 


\ip') = Itjj} , 


(94) 


which holds for any arbitrary state so that the norm of the eigenstate is given by 


(V’lV’) 


N,p 






(95) 


where we will choose the reference state \(j)) = |A^) 0 \ ® ■” ® order to keep the corre¬ 

spondence with the results obtained by Tschirhart and Faribault [18] transparent. By making use of the 
pseudo-deformation, the degeneracy of the deformed level is given by 2so(^) -I- 1. The normalization will 
be calculated for an arbitrary degeneracy, after which the limit of an infinite degeneracy follows from the 
contraction limit: 


ria a/ 1 + \/l + V'fi + o 

/i Q 1 li V ^ ^2 

{m = (</>! n = n + m{sT\ - soio), 

Oi a. 

{cl^m = (</*! 11*5^1^') = ^ det[J'^““'] (-so(C) + N\{Sf^o(^^-^\so{0) ■ 

u, 


(96) 

(97) 

(98) 


The matrices are given by Eq. (85), with the number of pseudo-deformed creation/annihilation operators 
No given by 2so(^), N and 2so(^) — N respectively. Since not only the expansion coefficients are necessary 
but also the matrix elements, these have been written out explicitly. 

Combining these expressions and evaluating the matrix elements, the result is 


N 

a=l 


(2so(C))!Af! det[J‘°*] 
{2so{0 - N)ldet[J^]’ 


(99) 


where the term (2so(C))!-^!/(2so('?) — N)l is obtained from the action of the generators on the irreps. Both 
and are m x m matrices with matrix elements 


jtot 

^ij 


2 ELi St +1 - St + +1)^ 


m-rij 


ifi = j 


( 100 ) 
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jdual _ 

— 


Vi-Va 


E N l-\-r)ir)a I 2 
Q;=l r)i — 

rii-Vj 


-E 




k^i,k^O rii-riu 


i^ + (2so(e)-iV + l)r7. if*=J 


3 


( 101 ) 


Taking the contraction limit of these expressions and making use of the prefactor introduced in Eq. (38), a 
determinant expression for the Dicke model is obtained as 


/,,,N .r, det[J*°*] 


with 


.n: 


det[J'^’' 


^j2E:=i^ + ^-Er^.^ if*=j 

3 


, ei-Cfc 


jdual _ 

'^ij ~ 


y^iV _1 


. gj-gP _ 1 if j = j 

ei — Xa ' G 2^k^i Ei—ek ^ ^ 


Hi 3^ j 


( 102 ) 


(103) 


(104) 


These correspond to the matrices derived previously [18] and can be written in the set of eigenvalue-based 
variables 

N 


A. = E^—. 

^ f:- — T- - 


(105) 


a=l 


eliminating any explicit dependency on the rapidities {a;^}- 

Again, a similar derivation can be made for the extended (p-I-ip)-wave pairing model (see Appendix B), 
resulting in 

(m = iV!SS ( 106 ) 


jtot _ 

— 


■det[J‘^““'] 

' ii+x„ _ Ho _ 1 V'" _ 1 _ n 

2- ^a= l ei-Xa, Ei 2 l^k^i ei-ek 2 ‘ ■> 

^ if , ^ j 


(107) 


and a dual matrix with an additional factor 1/2 in the first summation in the diagonal elements. Note the 
similarity between the diagonal elements of these matrices and the RG equations for both models. This 
similarity already arises at the level of the Gaudin algebra, where the RG equations are given by 


N 


1 + gSo{yZQa -f — Zia — g Z^a — 0, 


(108) 


i=l 


whereas the diagonal elements of the total matrix (J/°*) are (up to the term Z^i) given by 


N 


f 2 ^ ] Zia ~ 2so(0^iO ~ ^ ^ Zij — — 

^ Q=1 j/i ^ 


N 


1 + gso{C)Zoi + 2 ~ 2 ^ * ’ 

a=l 


(109) 


which resembles the RG equations, but with the roles of the energy levels and rapidities exchanged. 


6.3 Form factors 

Following the methods introduced in [17] and later used in [18] and [16], determinant expressions for several 
form factors of local spin operators can be straightforwardly obtained by making use of the results presented 
in the previous sections. These form factors for the Dicke model were previously given by [18], so we will 
restrict ourselves to form factors for the extended (p-l-ip)-wave pairing model. We only present final results, 
since the followed method is analogous to [16-18]. 
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6.3.1 Local raising and lowering operators 

The expectation values of local raising and lowering operators Sj, Si,b\b are non-zero between two states 
where the number of excitations differs by one. These can again be written by making use of the set of 
eigenvalue-based variables. For two states |{a:a}) and |{a;^}) with N and iV — 1 excitations respectively, we 
obtain 

({x„}|s;i{x„}) = ( 110 ) 

with an (rn — 1) x {m — 1) matrix with matrixelements 


f i gj+aic. I Y^(V-1 ej+Xf^ _ I 9„ _ 

k _ J 2 1 ei — Xa e^ —^ ei ' 




- 1 


-v/€^€fc 

ei-efe 


and defined as 


if f = j 
if i ^ j 


with i, j ^ k, (111) 


jdual _ 

^ij — 



+2^-E 


m €^+e 


Ci-ej 


i _ 1 


if i = j 

iii^j 


( 112 ) 


For the bosonic creation operator, we have 




with the m X m matrix given by 



ej+Xf^ _ „ 77 g I 9 _ Y^T? 

Z^u=l ei-x,, ^ a 2^.7 


Ci+e 


Oi-e 


i _ 1 


if I = j 
if f ^ j 


(113) 


(114) 


The expressions for the annihilation operators can be found as the Hermitian conjugates of these operators. 


6.3.2 Local ’number’ operators 

By making use of the Hellmann-Feynman theorem, the expectation value of local counting operators can 
be obtained. We will explicitly denote the eigenstates at a certain value of k as \xa{K)) and |a:^(«;)) with 
eigenvalue-based variables {A“} and {A^*} respectively. From the definition of the constant of motion Ri, 
the expectation value of can be obtained as 


{Xc,{K)\S°\Xa{K)) 



(115) 


with Ai defined by Eq. (54), where the set of these variables have to satisfy Eq. (56). By taking the partial 
derivative of these equations to k, a linear system can be found for as a function of A^, similar to the 
method originally introduced in [18]. Similarly, the off-diagonal expectation values can be found as 


(a;^(K)|S')!|a;„(K)) 



m 


Af)i: 

k=l 


dA“ detJ'= 
Ok det 


(116) 


with defined as 


= 


i(A“ + A(‘-E;^,f^-^ + 2«-l) ifz=j 


(117) 


Matrix elements for Wb can then be found by noting that N = b^b + 1^® expectation values 

of the counting operator are known. 
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7 Conclusions 


In the present paper we have shown how two different classes of integrable models can be obtained within 
one unifying framework, using the pseudo-deformation of a single quasispin. Starting from a trigonometric 
and a hyperbolic realization of the XXZ models, we were able to derive the integrability of the Dicke model 
and the extended (p + ip)-pairing model, respectively. Furthermore, this connection was then used to link 
the determinant expressions for overlaps, normalizations and form factors with the determinant expressions 
recently presented by us for XXZ RG models [16]. We were able to rederive the results for the Dicke 
model [18], circumventing the need for introducing the Algebraic Bethe Ansatz and the pseudovacuum, 
and generalize them toward the second class of integrable models containing a bosonic degree of freedom. 
Whereas strongly correlated systems in general exhibit unfavourable numerical scaling, expressions such as 
the ones presented in this paper allow for extensive numerical investigations of integrable systems. These 
results also further extend the description of RG integrable systems in terms of a new set of eigenvalue-based 
variables, explicitly symmetric in the rapidities, to models containing a bosonic degree of freedom. 
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A Equations for the eigenvalue-based variables 

For an XXZ model defined by a Gaudin algebra, it was previously shown [16] how the set of RG equations 
(6) are equivalent to a set of equations for the variables 

N 

A. = (118) 

a—l 

where we considered a system with n levels labelled i and N excitations labelled a. These equations are 
given by 

n ^ 71 N 

= riV(l - iV + 2^ - -A, + 2^ (A, - A,) + (1 - 2^,) ^ (119) 

j^i ^ j^i a—1 

where F is a constant defined by the Gaudin algebra as 

Xl - 4 = F, Vz ^ p (120) 

It can be seen that these equations are closed in the set of variables {Ai} if Si = Starting from an 

{m-\- l)-level system and introducing the pseudo-deformation of a single level 0 and taking Si = ^ 0, 

a closed set of equations can similarly be obtained. The equations for the set of A^ are then given by 

A2=FV(m-iV + 2^)--A,+ ^ (A, - A,) + 2^Z,o(A, - Aq). (121) 

By introducing the proposed parametrizations for both models models, a closed set of equations can be 
obtained in the contraction limit. 


A.l Dicke model 

Starting from the trigonometric representation (F = 1) for the Dicke model, we have 

N N 


Ao = 

Ol — 1 


a. = E 


1 + Pi'na 
Vi - Va 


i = 1... TO, 


( 122 ) 
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which have to satisfy 


N 


1 + rjiTja 

.h 


o 1 4- r 

=Nim -N + 2so/0 --Y1 


a=l 


ViVa 
Vi - Va 


39^i 


1 + ViVj 

V^ - Vj 


' N 

E 


1 + ViVa 1 + VjVo 


N 


Vi - Va Vj - Va 


„S0 / l+V^Va 

^ Va=l 


(123) 


Introducing the proposed expressions for r]i and g, multiplying everything with 2so/^, and taking the 
contraction limit immediately simplifies this to a set of equations 


™ 

=N- Af (e. - eo) + E 


(124) 


in the variables 


A: 


D 


N 


E 



(125) 


A.2 (p + ip)-wave pairing 

Starting from the hyperbolic parametrization (F = —1), we have 


Ao - E 


+ 2soXo 
^Vo - 2soXo 


A = E 


Si + Xo 
Si - Xo 


(126) 


where the following expansion holds for ^ 1 


N 


Ao = E (-1-^^) + 0(^2) 

\ -Qr. T- / (in 


a—l 


Sq Xq 


C ^ ^ 

^ ^ Vo 


So Xo 


(127) 


where we can define Aq = X]c(=i Vo/^a- Introducing this parametrization in the set of equations, we obtain 


A- = -N{m-N+‘^)-2 ( y + «:) A^ + E 


So 




+ Cj 

e^ - e. 


2so f SsoCi + ^77j 




(A,-A,) + — (A, + Ar+_Ag)+0(e), 

(128) 

where all terms containing 1/^ drop out. In the contraction limit, this results in 

A? = -Nim -N)- 2 kA, + E - A,) + ^(A + N) + 2Ag. (129) 

An additional equation can be found linking Ag to the set of A^ by considering (119) for the pseudo-deformed 
level. Keeping only the terms to dominant order in we obtain 

= -Nil-N+m)-2 i-N- 1-Ag)+E(lV+A,)+ (^1 - (^iV + ^ + Ag^ +0(0, (130) 

where all terms containing 1/^ again drop out. In the contraction limit we are left with 


-2AS + y: Aj = 2 kN, 
1=1 


( 131 ) 


where we have obtained a closed set of equations in the eigenvalue-based variables. 
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B Results for the {p + ip)-wa.ve pairing model 

In this Appendix, it is shown how overlaps for the extended (p + ip)-wave pairing model can be derived in 
the contraction limit of a hyperbolic RG model. The Bethe Ansatz state is given by 



= E n (^ 1 ) ’ 1 ^) (132) 

l{N^}] i=l 


for the parametrization proposed in Section 4.2. The expansion coefficient for a state |fVo; {*a}) is given by 


lim TVT 

C-fO TVq! 

(¥: 

«0 

per 

/ ^Oqi 
-^0Q2 

■ ■ ■ ^^0CK2 

^l\Oi2 

^IN-NqCii \ 

^IN-NqO/.^ 




^ -^OCKiV 

■ ■ ■ -^Oaiv 




No 


(133) 


with Xij = 2yAi^/(ei — ej) and eo = ^r?g/(2so)- Introducing a different deformation parameter for each 
column, and relabelling eg = ^/( 2 so? 7 o) Ci i^ column i, this permanent can be rewritten as 


^[{N^}] = -]vT 




lim ... lim det J, 


Vo°Nol ynEi^“ Ca 

with the N X N matrix J again consisting of four blocks, defined as 


with 


Jij 

J^a 

Jai 

Jab ^ 


Jij — 


Jab — 


J = 


' Y^N Q+Xc, _ y^No Ci+ik _ s^N-No Ci+gg _ 1 ;f 

A^a=l l^k^i l^a=l Ci-Eg u i — _y 

2aE 


'C,-0 


E N gg+Xg iVO eg+gg gg + t;^ 1 '£ _ 7 

OL — 1 ea — Xoc 2-^C^a gg —gg 2-^i — l Ea—Qi 

2,/g^ 


^N — Nq ga+€c 


''-^0 gg + G 


eg-Cb 


iii^j 

if a = 6 
ii a ^ b 


Jia — Jai — 


V Ci^a 


G-ea’ 

where we have used the notation = £„. Taking the limits to zero as 


Cl < C 2 < ■ • • < Cwo: 


(134) 

(135) 

(136) 

(137) 

(138) 

(139) 
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this results in 


J,, = -iV - ^ ^ + {N-No)-l 


Ci + 0 


i<i 




= -iV - (i - 1) + {No -i) + {N - No)-l = -2i 


Jrj =0 


N N—Nq 

j 

Jaa — >->--/Vo - 

^ ea- Xc ^ Cn - Ch 


Jab — 


a=l 


b^a 


ea - ^b 
Jia — Jai — 0- 
The expansion coefficient is then given by 


(140) 

(141) 

(142) 

(143) 

(144) 
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which implies that the overlap of the wavefunction 
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with a non-interacting state \(j)) = |A^oi {*a}) is given by 
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The expression for the normalization can immediately be taken from the Dicke model, so we have 

The off-diagonal elements of both matrices are given by l{^i — while the diagonal elements of the 
top matrix are given by 
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and the diagonal elements of the bottom matrix are similarly found as 
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